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Abstract 

It is well-known for vibro-impact systems that the existence 
of a periodic solution with a low- velocity impact (so-called 
grazing) may yield complex behavior of the solutions. In 
this paper we show that unstable periodic motions which 
pass near the delimiter without touching it may give birth 
to chaotic behavior of nearby solutions. We demonstrate 
that the number of impacts over a period of forcing varies 
in a small neighborhood of such periodic motions. This al- 
lows us to use the technique of symbolic dynamics. It is 
shown that chaos may be observed in a two-sided neighbor- 
hood of grazing and this bifurcation manifests at least two 
distinct ways to a complex behavior. In the second part 
of the paper we study the robustness of this phenomenon. 
Particularly, we show that the same effect can be observed 
in "soft" models of impacts. 

©2012 L&H Scientific Publishing, LLC. All rights reserved. 



1. Introduction 

Vibro-impact systems appear in different mechanical problems (modeling of impact dampers, 
clock mechanisms, immersion of constructions, etc.). Their properties show a lot of resem- 
blence to classical nonlinear systems. Particularly, chaotic dynamics is possible [1-21]. 

The so-called grazing bifurcation was first described by A. Nordmark [6l [17] . The critical 
value of the parameter corresponds to a zero velocity impact of a periodic solution. It was 
demonstrated that this bifurcation implies non-smooth behavior of solutions, instability of 
the periodic solution in the parametric neighborhood of grazing and, in additional assump- 
tions, chaotic dynamics. It was shown theoretically, how periodic solutions corresponding 
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to different values of impacts (ni, n2, . . . , rim) over successive periods of forcing may appear 
p]. 

However, parametric neighborhoods of the so-called continuous grazing [T71 118] consist 
of two parts. One of them corresponds to periodic solutions with a low- velocity impact, 
another one corresponds to periodic motions that pass close to the delimiter but do not 
touch it. In our paper we consider the second case assuming that the mentioned periodic 
motions are unstable. 

We consider the system studied in [H]: the variable x G [—1, 1] satisfies the equation 

x + 2/3x-x = fit) (1) 

over free flight intervals. Impacts, given by the equalities 

f + = — rf _ where v± = x{t ± 0) 

take place at |x| = 1. It was assumed that the right hand side of the equation can be 
represented combination of Dirac functions, 

oo 

/W=7 Yl {6{t-kT-a)-6{t-kT)), 

k=—oo 

where a G (0,T). For the considered system the authors proved the existence of so-called 
non-classical bifurcations, corresponding to coincidence of the impulse action which takes 
place for t = a + kT and t = kT and the impact. These bifurcations are accompanied 
by the disappearance of the periodic solution and, sometimes, by appearance of a "strange 
attractor". 

The main aim of our article is to show that chaotic dynamics is possible in a two-sided 
neighborhood of grazing. Namely, it can be caused by the presence of a unstable periodic 
motion passing near the delimiter. An experimental illustration of this is given in [T6] . 

We start by considering an impulse model of impacting system where interaction between 
the moving particle and the delimiter is instantaneous. However, from a physical point of 
view the so-called "soft" model [21 [TOl [TTl [16] where the delimiter is considered as a very 
stiff spring, is more relevant. We prove that any hyperbolic invariant set that appears in 
the impulse model persists if we pass to a "soft model" with a sufficiently stiff spring. 

The rest of the paper is organized as follows. In Section 2 we give a description of the 
impulse model of impacts. In Section 3 the grazing family of periodic solutions is considered. 
Section 4 is devoted to description of a set of initial conditions corresponding to solutions 
with a zero velocity impact. In Section 5 the main result on the existence of chaotic invariant 
set is formulated. In Section 6 asymptotic estimates for Jacobian matrices corresponding to 
near-grazing solutions are given. The existence of a transverse homoclinic point is proved 
at Secition 7, and in the next section the corresponding symbolic dynamics is described. 
A simple example of a single degree-of- freedom system is considered in Section 9. Also we 
prove in Section 10 that an infinitely stiff delimiter may be replaced with a sufficiently stiff 
spring. The conclusion is given in Section 11. 
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Later on we use the following formalism: indices of parameters of successive impacts 
of a fixed solution (phases, velocities, etc.) are denoted by superscripts in order to distin- 
guish them from ones of coordinates of vectors always denoted by subscripts. Denote by 
col(ai, . . . , am) the column vector, consisting of elements ai, . . . , am- Here a^ themselves 
may be vectors. Also, we use the notation a = col(a2, . . . a„). 

2. Impulse model 

We study the motion of a point mass, described by system of second order differential 
equations of the general form and impact conditions of impulse type. 
Consider a segment J = and a smooth function 

f{t,z,fi) :R2"+^ X J^M". 

Here is a scalar parameter. For example, this can be the amplitude of free periodic 
oscillations of the considered system or the coefficient of restitution. Suppose that 

f{t,z,ij,) = f{t + T,z,fi). 

Let 

z = co\{zi, . . . ,Zn); Zk = col{xk,yk), k = l,...,n; z = co\{z2, . . . , Zn); 
/ = col(/i,...,/„); a; = col(a;i,...,a;„); ?/ = col(?/i, . . . , 
X = col(a;2, . . .,Xn), y = col(?/2, • • •,!/«)• 

Consider the system 

Xk=yk] yk = fk(t,z,fi), k = l,...,n. (2) 

Let the coefficient of restitution r = r(yu) G (0, 1] be a - smooth function. Suppose 
that Eq. (|2]) is defined for 2; G A = [0, +00) x and the following impulse type impact 

conditions take place as the component Xi vanishes. 

Condition 1. 

1. If xi(to) = then x{to + 0) = x{to - 0), 

yi{to + 0) = -r{f^)yi{to-0), y{to + 0) = y{to - 0) . (3) 

2. Xi{t) ^ for all t where z{t) is well-defined. 
Consider the vibro-impact system 

i = y; y = f{t,z,n); 

Condition 1 is applied if xi = 0. 
We give two definitions of solutions of vibro-impact systems. 

Definition 1. The function z(t) = co\{xi(t),yi{t), . . . ,Xn{t),ynit)) is a solution of Eq. 
(j4]) with a finite number of impacts over the interval (a, b), if there exists a finite number of 
instants a = to < ti < . . . < tN < tN+i = b such that the following conditions are satisfied. 
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1. All components of the solution z{t), except yi{t), are continuous while t G {a,b). The 
discontinuity set of the function yi(t) is a subset of {tfc}. 

2. The function Xi(t) is non-negative on (a, b). The set {ti, . . . , t^} is the set of all zeros 
for this function. 

3. For any k = 1, . . . , N 

yiitk + 0) = -r(/i) yiitk-0). 

4. The function z{t) is a solution of (|2]) on every interval 

For completeness of the mathematical model, we define a solution with an infinite number 
of impacts. 

Definition 2. The function z{t) = col(xi(t), ?/i(t), . . . ,Xn(t),ynit)) is a solution of the 
vibro-impact system (jl]) over the interval (a, b) if there exists a disjoint sphtting (a, 6) = 
1+ [jXo [JX_ with the following properties. 

1. The set X+ = {t G (a, 6) : > 0} corresponding to free flight motions is open. The 
set I_ = 

{te{a,b):xi{t)=yi{t) = 0,fi{t,0,0,X2{t),y2{t),...,Xn{t),yn{t),fi)^0} 
is closed. 

2. All the components of the solution z{t), except yi{t), are continuous over {a,b), the 
discontinuity set of the function yi{t) is a subset of Iq. 

3. The function Xi{t) is non-negative on (a, b). The set Xq IJ-^- zeros of this 
function. 

4. The function z{t) is a solution of (|2]) on every open convex subset of X+. 

5. The set Iq is finite or countable. All limit points of this set belong to X_. 

6. For any to ^ -^o Eq. ([3]) holds true. 

7. The function z(t) is a solution of 

Xk=yk] yk = fk(t,0,z,fi), k = 2,...,n 
on every open convex subset of X_. 

Generally speaking, solutions with infinitely many impacts are not unique. For example, 
one may take the single degree-of-freedom equation x = —1 with delimiter at a; = and 
r < 1. The solution x = is non-unique "backwards". 

We identify the vibro-impact system with the pair (/, r). 
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Introduce the topology on the set X = X{J, n, T) of vibro-impact systems, corresponding 
to fixed J, n and T. This is the minimal topology such that for all (/o, ro) € X{J, n, T) and 
every R> the set 



+ 



{{f,r) e X : sup(i_^_^) {\f{t,z,fi) - /o(t,z,/i)| + 



df , . dfo, , 
|r(/i) - ro(/i)| + 



(9/ , , (9/o , - 
r'ifi)-r'M\)<R} 



+ 



is open. The space A* with this topology is Hausdorff. 
3. The family of periodic solutions 

Since the solutions of the vibro-impact systems are discontinuous at the impact instants, 
the classical results on integral continuity are not applicable. Nevertheless, the following 
statement is true. 

Lemma 1. Let z{t) = co\{xi{t),yi(t), . . . ,Xn(t),yn(t)) be the solution of @ for fi = fi* and 
initial conditions z(tQ) = Zq = col(a;jQ, yj'o' • • • 5 ^no? 1/no) ■^^c^ ^^'^^ ^10 0- Suppose that 
this solution is defined on the segment 3 to- Assume that there are exactly N zeros 

t^ < < . . . < t'^ < t+ of the function xl{t) over the segment [t^,t^] and yi{T^ — 0) 7^ 0, 
(j = 1, . . . ,N). Then for any e > there exists a neighborhood U of the point co\{zq, ^*) 

N 

such that for any fixed t G \ |J {t^ — e, r'^ + e) the mapping z{t,to, zq, ft) is - 

fe=i 

smooth with respect to the variables (tc-Zo,/^) ^ U. Moreover, these solutions have exactly 
N impact instants T^{to, z^, /i) (j = 1, . . . , A^) over the segment t+]. These instants and 
corresponding velocities 

-yiiT^to, Zq, Hi) - 0,ti,z\/ii) 
- smoothly depend on tQ,ZQ and fiQ. 

Proof. Let the number k be such that to ^ [t^,t^^^] (assume, if necessary r° = t_, 
T^'^^ = t^). The solution z{t) of (jl]) is also one of (|2]) over any free fiight segment. The 
impact instants r*^ and t''^^ as well as the impact points for solutions z{t'' — 0,to,ZQ,fi) 
and z(r'^+^ — 0,tQ,ZQ,fi) smoothly depend on their parameters. Similarly, values r*^"^ and 
z{t^^^ — 0,t(),zo,fi) smoothly depend on and z{t'' — 0,tQ, zq, fi), as well as r'^"*"^ and 
2^j-k+2 _ o,to,-2o5yu) are - smooth functions of r^^^ and z{t''^^ — 0,to,Zo,fi) and so on. 
□ 

Suppose that for /i ^ /^o the considered system has a family of periodic motions, that 
pass near the delimiter, and touch it if and only if = /io (Fig. 1). More precisely, the 
following condition is satisfied. 

Condition 2. There exist a segment J = [yUo,/Ui] C J and a family of T - periodic 
solutions 

ip{t,fi) = C0\{(pxl(t,fl),(Pyi{t,fl),...,(p:^n{t,fi),(Pyn{t,fl)), t E R, fiEJ 

of (jlj) with the following properties. 
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Fig. 1: The grazing family of periodic solutions. 

1. For any pair {to, ft) G [0,T) x J, such that {pxiito,fL) ^ the function (p{t,fi) is 
continuous in a neighborhood of the point {to,fi). 

2. For any e J the component ipxi{t,fi) has N zeros t^{^j) < . . . < T^(/i) over the 
period [0,T). The values r'^(yu) and the velocities 

Y'if^) = -^y^iT\fx)-0,fx), k = l,...,N 

continuously depend on fi. 

3. For fi > fiQ the component {pxi(t,fi) does not have any other zeros. The function 
^xi{t,^io) has exactly one additional zero t = 0. 

4. Velocities Y''{fi) are such that 

'fyiiO, yWo) = 0, /i(0, 0, 0, (^(0, /io), /io) = 0o > 0, 

F^(/i)>0, V/iGJ, fc = l,...,iV, (5) 
min r^(/i) > 0, max r^(/i) < T. 

Here ^{t, ^) = C0\{(px2{t, fl),(Py2(t, fl), . . . ,(p^n{t, fi),iPyn(t, fi)). 

If /i is the amplitude of free periodic oscillations of the system, the critical value fiQ 
corresponds to the distance between the neutral position of the particle and the delimiter. 
However, in this case, to fulfill Condition 2 we must replace /i with — /i. 

Later on we may suppose without loss of generality that <^(0,/io) = 0, /io = 0. 

Denote 

Vx{t, IJ) = C0l(<^a;i(t, /i), . . . , ^xn{t, /i)), 
iPy{t,^) = C0\{(pyi{t,^l),. . .,ipyn{t,n)), 
(Pxit,fl) = C0\{ip^2{t,fl), ■ ■ .,<fxn{t,IJ')), 
(fy{t, fJ,) = C0l{ipy2{t, /i), . . . , iPyn{t, fJ,)). 
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Consider the shift mapping for system (jlj) given by the formula 

Sf,{zo) = z{T + 6, -9, zq, /i). 

Here the value 6 = 6{fi) will be specified later (proof of Lemma Hj). 

For small positive n the mapping is - smooth in a neighborhood of its fixed point 

4. Separatrix 

Denote = {^o e A : 3ti G [-T,T] : 0, ^o, = 0}. 

Lemma 2. There exists a neighborhood Uq of zero such that if the parameter ^ is small 
enough, the set f/o zs a 2n — 1 dimensional surface, that is the graph of the smooth 
function xi = -y^{x,y) (Fig. 2). Moreover, 

where 7^ zs a smooth function such that 7^(0) = 0. 




Fig. 2: The near-grazing behavior of solutions. 

Proof. Take a point po ^ T/i- Let to be such that 2:1(^0, 0,po) = 0, s = t — to, 

2;i(t + 0,0,po) = col(xi(t),yi(t)). 

Let us show that if to is close enough to 0, we may take so ^ |to| so that the function 
Xiito + s) does not have any zeros on [— so,So], except s = 0. Otherwise, there exists a 
sequence t° ^ (suppose without loss of generality, that t° > and the sequence decreases), 
a sequence t\ G [0,t°) and one, consisting of solutions: 

COl(Zfci(t), . . .,Zkn{t)) = C0l(Xfcl(t),yfei(t), . . . ,Xkn{t),ykn{t)) 
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Fig. 3: Sequences t^. 

of the system (jlj), such that Zkiit^) = 0, Xkiitl) = (Fig. 3). 

Also, there exist time instants tl G (t^, t°) such that iki{tk) — instants tl G (tl, t^) 
such that Xki{t\) = 0. Moreover, t| — )■ 0, Xki{tl) — )■ 0, Xfci(4) — ^ 0. Without loss of 
generality, we assume that z{tl) — )■ zq. Then — )■ /i(0,0,/i) = 0. This contradicts (ISl). 
Then for all s G [— sq, sq] the function xi(to + s) can be represented as series 

xi (to + s)= X2S^ + Xss^ + . . . (7) 

Differentiating (I?!), we obtain that xiito + s) = 2X2S + 3X35^ + . . .. On the other hand, 
X2 = Xi(to + 0)/2^/i(0,z^,/i) as to^O. Then 

Xi(0) = /i(0,2:^,/x)(to)(l + o(l))/2; 
2/1(0) = /i(0,z^,/x)(to)(l + 0(1)). 

Since yi = Xi, is true. □ 

Take a small parameter <^ > such that the sets 

V^ = {zeA: \\z-zj ^ c^} C f/o, 

V~ = {{z eV^:xi< -f^{x,y)},V^ = {{z e V : Xi > 7^(^,1/)} 
are correctly defined and nonempty. Here || ■ || is the Euclidean norm. 

5. The main result 

Consider the (2n) x (2n) matrix 

A = (aij) = (Ai, . . . , A2n) = lim 05^(2^). 

Suppose the following condition is satisfied. 

Condition 4. The matrix A does not have any eigenvalues on the unit circle. At least 
one of the eigenvalues of A is inside the unit ball and at least one is outside this ball. 

Note that the matrix A corresponds to the motion out of a neighborhood of grazing. 
This may be a motion, described by a linear system without any impacts. 
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Fig. 4: Homoclinic points. 

Let be the linear hull of eigenvectors, corresponding to eigenvalues inside the unit ball 
and M" be the space, corresponding to eigenvalues outside the unit ball. Let k = dimM*, 
TTi be the hyperplane of the delimiter, defined by the equality xi = 0. 

Condition 5. M^M" ^ tti. 

Then the intersections of these spaces with a hyperplane tti are transverse. Denote 
TT*'" = TTi f]M^''^. Let aij be entries of the matrix A^'^. Select a basis el, e^, . . . , in the 
space and one e", 65, ... , e2„_fc in the space M" so that eJXvrJ, e"±7r"; e| G vrj, e" G vr" 
for j > 1. Denote the components of vectors ej by e^j, a G {s, -u}, i, j = 1, . . . , 2n. Both the 
values ell nonzero. Denote 

i?^ = e^iMi, i?" = e^i/e^'i, i?'^ = a22/ai2, i?" = «22/«i2. 
Condition 6. Either 

ai2 > 0, (i?'' - R')/{R' - W) > 0, (8) 

or 

ai2 > 0, (i?" - - i?") > 0. (9) 

From the geometric point of view. Condition 4 means that the point 2;^ has nontrivial 
stable and unstable manifolds {VV^ and respectively). Condition 5 means that for small fi 
these manifolds intersect transversally the surface and bend at the points of intersection. 
Condition 6 means that the "prolongation" of the manifold beyond the intersection with 
F^ intersects transversally with W or vice versa (Fig. 4). 

Later on we suppose that ([8]) takes place. Otherwise, if ([9]) is true, we consider the 
mapping instead of 5^. In the proof the matrix A and the corresponding elements 
should be replaced with the matrix A~^ and corresponding elements and all the references 
to ([8]) should be replaced with ones to ([9]). 

Theorem 3. Let Conditions 2 and 4 — 6 be satisfied. Then there exists a value Hq > 0, such 
that for all fi G (0, fio) the mapping has a Devaney chaotic invariant set f2^ . 
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Namely, there is a hyperbolic transitive infinite invariant set of the mapping S*^ such 
that periodic points of S*^ are dense in Kf^. 

6. Grazing 

Now we start proving Theorem [3l Fix a small value /i > and a solution 

z*it) = co\{xlit),ylit),...,x:{t),y:{t)) 

of the corresponding system, having an impact at t = to- Suppose that the corresponding 
normal velocity Iqi = ~l/i(^o ~ 0) is nonzero but small. Denote Zq = z*{tQ — 0). Fix a 
positive value Sq and consider the mapping G{qo) = z{to + So,tQ — so,qo,fi), defined in a 
neighborhood of the point = -2* (to — sq). Here we assume that the point and the 
parameter sq are chosen so that there exists a neighborhood Q 3 qQ such that any solution 
z{t) = z{t, to — Sq, z_, fi) (z_ G Q) has exactly one impact over the segment [to — Sq, to + Sq]. 
Denote the corresponding instant by ti = ti(-Z-) and the normal velocity, defined similarly 
to Foi, by Yi = Yi{z-). Let Z = z(ti) be the tangent component of the solution z(t) at 
the impact instant. Take the values s± = s±{z-) so that to ± sq = ti{z^) ± s±{z-) for all 
z_ G Q. The mapping G is smooth in the neighborhood of the point qQ, let us estimate the 
Jacobian matrix DG{qQ). Let 

Z- = C0l(a;i yi . . . , Xn-, Vn,-), 
"JC — col dC ^ ^ • • • ^ OC ^ J 

y_ = col(yi _, . . .,?/„-), 

X_ = C01(X2_, . . . 

y- = co\{y2,-,...,yn,-) 

Denote 

z+ = z{to + So) = z{to + So, ti + 0, 0, rYi, Z, /i), 
x+ = x(to + So) = x{to + So, ti + 0, 0, rVi, Z, fx), 
y+ = y{to + So) = y{to + so, ti + 0, 0, rYi, Z, /i). 

Similarly, we define the values Xi^+, yi,+, x^, y+. Consider the Taylor formula for values 
z± as functions of s±: 

a;i,_ = Yis_ + f,{ti, 0, -Fi, Z, fi)sl/2 + p,,4s^,ti, Yi, Z, i2)sl; 
yi- = -Yi - /i(ti,0^-Fi,Z,/i)s_ + Py_(s_,ti,yi,Z,/i)s_;_ 
x_ = x(ti) - y(ti - 0)s_ + /(ti, 0, -Fi, Z, fi)s^_/2 + p5,-(s_, ti, Yu Z, fi)sl; 

y- = yiti -0) - /(ti,0, -Fi^Z,/i)s_ + p5_(s_,ti,Fi,Z,/i)s_; , . 

= rY^s+ + /i(ti, 0, rFi, Z, p)4/2 + p,,+ (s+, ti, Fi, Z, ^ ^ 

yi,+ = rYi + /i(ti, 0,_rYi, Z, /i)s+ + py,+ is+, ti, Yi, Z, p)s+; _ 
x+ = x{ti) + y{ti - 0)s+_+ /(ti, 0, rYi, Z, p)sl/2 + ps,+ (s+, tj, Y^ Z, /x)s^; 
y+ = yiti - 0) + /(ti, 0, rYi, ^, p)s+ + Pg,+ (s+, ti, Fi, Z, p)s+. 
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Here all functions, denoted by the letter p, are smooth with respect to all arguments 
except s± and vanish as s± = 0. Denote 

/ofc+ = fk{to, 0, rloi, Zo, fi), fok- = fk{tQ, 0, —Yqi, Zq, p). 

It follows from (fTOl) that 



dz4 







\ 




/oi+ 


r 


s+=0,Yi=Yqi,Z=Zo 


W+ 


-^271-2 / 




( >01 









-foi- 


-1 


_=O^Yi=YouZ=Zo 


^ Q- 


-E'2„_2 



d{s^,Y,,Z) 

Here i?2n-2 is the unit matrix of the corresponding size, 

Q+ = C01(?/2(^0 - 0), /o2+, • • • , Vnito ' 0), fon+), 



Q- = Col(-|/2(to - 0), -/02+, • • • , -Vnito - 0), -/o„+). 

Denote 

r, _ dfk{to,0,yi,Zo,n) 

Jk — 



dyi 



yi=o 



Clearly, ds^/ds^ = — 1. Then, similar to the results of the paper [TO], we obtain 



B = lim — — ! 

s±^o dz- 



— r 








hi 


— r 





Bi 








(Fig. 2). Here 5i = col(63i, • • • &2ni), 

b2i = -(/oi+ + r/oi-)/roi = -(r + l)0o(l + O(roi))/>^oi, 
fc2i-ii = 0, 62,1 = {foj+ - foj-)/Yoi = (r + 1)/^ + 0{Yoi). Note that det B = r\ 
7. Homoclinic point 

Clearly, the mapping is smooth at the points of the set V^, except ones of the curve F^. 
The matrix DS^^zq) is of the form A + -R+(zo5 /^) where 

sup -R+(2:o5 /^) ~^ 

as p — )■ 0+. For zq G there exist values to(-2o) and Yqi = ^oi(-2^o), continuously depending 
on 2:0 and such that 



xi(to,0,2;o,/i) = 0, t/i (to - 0,0,^0,/^) 



-^01, 
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where loi(zo) — ?■ as dist(zo, T^) — )■ and 




The matrix D = 05^(^0) can be represented in the form yl^(zo)-B^(-2o)- Here 



A^{zo) = A + R_{zo,fi) 



sup R-{Zq, /i) — 7- 



as /i — )■ 0+. The matrix B^{zq) can be represented in the form 

B^{zo) = {E + Rb{zq, ^i))Bq^^{zq), 

where RB^ZQ.ji) = 0(1) as ;U — )■ 0+ uniformly by Zq G V~ and the matrix Bq^^j^^Zq) has 
the form (ITT]) . Suppose A' = A^(2;o)(-E' + Rb{zq,h)). Denote columns of this matrix by A'- 
(j = 1, . . . , 2n). Then 



Since deti5(2;o) = detA^i^zo) detB^^zo) = (r^ + 0{^))Aq and ai2 7^ 0, one of the eigen- 
values of the matrix D{zq) is equal to 



The corresponding eigenvector is of the form m+ = A2 + 0{^). The linear space, correspond- 
ing to other eigenvalues, tends to the hyperplane tti = {(0, 1/1,2)} as /i — )■ 0. Since 012 7^ 0, 
the vector u+{zo) does not belong to vri. The vector A2 (as well as the vector A'2) is out of 
the space and the space M". 

The matrix D~^{zo) satisfies the following asymptotic estimates: 



D{zo) = A^{zo)B^{zo) = 

(-(r + l)A^0o(l + O(Foi))/Voi, 

-rA',{l + O(Foi)), ^3 + O(yoi), • • • , + OiYo^)). 



{r + l)auMi + 0{Yoi))/Yoi. 



I 



-rA-^ +0(/i) 

(r + l)Ar(l + 0(/x))/roi 
A^ + 0{fi) 



\ 



D 



1 



v 



^2n + ^(/i) 




of the matrix D ^. The corresponding eigenvector is m_ = 62 + 0{fi). Here 



62 = col(0,l,0,...,0). 
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It follows from the Perron theorem that in a small neighborhood of the point Z/j, there ex- 
ist local stable and unstable manifolds Wf^^ and W^^ respectively. Both of them are smooth 
in a neighborhood of the point z^. Clearly, M*'"" = T^^Wi^^. Extend the manifolds Wf^^ 
and Wl^^ up to invariant sets of the mapping S^. Denote these sets by and W^, respec- 
tively. Generally speaking, both of these sets consist of a countable number of connected 
components. All these components are piecewise smooth manifolds. 

Lemma 4. The intersection of the sets and contains a point p 7^ (Fig- 4)- The 
connected components of intersections of the sets and VF" with a neighborhood of the 
point p, containing this point, are smooth manifolds and their intersection is transverse. 

Proof. If the parameter ji is small enough, the manifold VF" transversally intersects 
with the surface P| V^. Denote the k — 1- dimensional smooth manifold, obtained in the 
intersection, by g". 

Let = S^{q'^). The neighborhood of the point may be chosen so that both 
the manifolds and g" intersect with V^. Let L" be the connected component of the 
intersection of (W^" P| V^) \ (g" |J g") whose boundary intersects with g" and g". We select 
the parameter 9{ijl) in the definition of so small that g" and L" are correctly defined. 

The neighborhood may be chosen so that diam ^ as /x — ?> 0+. Let = S^{L^). 
For any z E L'^ the tangent space M'^{z) = T^W^ is the linear hull of unit orthogonal vectors 
Vi{z),V2{z), . . . , V2n-k{^), taken so that for all /c = 1, . . . , 

lim max \\v'^(z) - el\\ = 0. 

Recall that vectors have been defined immediately after Condition 5. 

The surface is not smooth in a neighborhood of the manifold g". For the points z G 
L" the tangent space M'^{z) — T^W"^ is a linear hull of unit vectors v'({z),V2{z), . . . , v^^{z) 
such that 

lim max lli'tf^;) — 6^11 = 

for all A; > 1 and 

lim niax||i5]'(z)->l2/||^|| 11=0. 

Here we use the fact that contains the manifold g" whose inclination with respect to tti 
is arbitrarily small if is small. 

It follows from Condition 6 that for all z E L\ vectors 

vl{z),vl{zl...,vlSz) 

are linearly independent. Vectors A2 and e" belong to different half-spaces, separated by 
the hypcrplane vri. Consequently, for all zq G L", Zi G and Z2 G W^" H ^ vectors Vi{zq) 
and (2:1) belong to different half-spaces, separated by the linear hull of vectors 

^2(^2), . . . , VI^{Z2),V{{Z2), . . . , <(^2). 
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Let Pfj^ be the nearest to point of the set tti f] M". The second coordinate of equals 

to 

yo = yi,^ - (12) 

Recall that Xi^^ and yi^^ are the first two coordinates of the point z^. 

Condition (|8]) is equivalent to the existence of a solution (X, Y) = {xp, i/p) of the system 

Y - ' Y -yo 

such that Xp > Xi^^. Here the value z/q may be found by (fT2l) . Hence, the affine space, 
containing the point po a^^d parallel to the linear hull of the vectors A2, 63, . . ., intersects 
the space, tangent to the manifold W"^ at the point z^. The first coordinate of the intersection 
point is greater than xi^^, that provides a transverse intersection of L" and W^, if /i is small 
enough. □ 

The parameter /i is chosen so that the point of this intersection, nearest to 2;^ does not 
belong to the hyperplane tti. 

The Smale-Birkhof theorem [23] on the existence of a chaotic invariant set in a neigh- 
borhood of a homoclinic point is not applicable in the considered case since the mapping S"^ 
is discontinuous. So we have to find a chaotic set of the mapping S*^ "manually". 

8. Symbolic dynamics 

Consider new smooth coordinates col(Ci, . . . , Cn^; C") • • • ; Cn„) the neighborhood of the 
point z^ such that the following conditions are satisfied. 

1. The point z^ corresponds to zero in the new coordinate system. 

2. The Euclidean norm of any column of the matrix 

gcoi(c^...,c,cr,---,cj , , 
Fz ^'^^ 

equals to 1. 

3. In a small neighborhood of zero the stable and the unstable manifolds are given by 
the conditions Q = and Q = {j = 1, . . . ,n) respectively. 

4. The direction of the tangent line to the axis, 0(1, coincides with one of the vector 
and one, corresponding to the axis OQ coincides with the vector U-. 

Consider the neighborhood Qq of the point 2;^, defined by conditions 

iqi^ss, igi^su, j = i,...,n 

(Fig. 5). For any m G Z define Qm = S^{Qo). Denote parts of the boundary of the set Qq, 
that correspond to Ci = ±£^s and C" = by and dy respectively. Here the positive 
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Fig. 5: Domains Hq and Hi. 

values Eg and e„ are chosen so that Ss ^ and there exist natural numbers m+ and m_ 
such that p G Qm+ flQ-m- (Fig. 6). We may take n and, respectively, eigenvalues X^{Zfj) 
so that C V~ for any — m~ < j < m~^. 

Due to Lemma m the set Q-m~ r\Qm+ contains at least two connected components. One 
of them denoted by Hq contains the point 2;^, another one denoted by Hi contains the point 
p. Let m = m++m-,Ho = S^-{Ho), Hi = S^-{Hi), H = H^\jHi (Fig. 5). 

Let us show that the set 

= fl S^\H) 

k€Z 

is chaotic. Evidently, the set K is invariant with respect to the mapping S^, compact and 
nonempty, since it contains the point 2^. Moreover, neither the inverse images S'^^tti 
{k G Z) of the hyperplane tti nor ones of the surface F^ intersect with K. Consequently, for 
any integer k there is a neighborhood f/^ of the set K such that the mapping S^lu,, is 
smooth. 

Lemma 5. For any A; G N, any set i = (io, ■ ■ ■ ,ik) such that ij G {0, 1} for any j = 0, . . . , k, 
the set Hi = Hi^ f] S^^^^Hi,) H • . • fl S-'^^Hi^) is not empty. 

Proof. Consider an arbitrary arc r], joining disks and d~ and defined as the graph 
of a smooth function 

(Cl; • • • ; Cn5 C2 5 ■ ■ ■ ? Cn) — ^(Cl ); 

such that max [[^'(Ci*)!! ^ 1- Let us call such arcs admissible. Similarly to the Palis lemma 
[23l chapter 2, Lemma 7.1], one may show that for small values of /x. Eg and e„ there exists 
an embedding of the curve S~"^{r]) to the manifold W^, arbitrarily C^-close to the identical 
embedding of S~"^{ri) to M^". Particularly, this means that the set S~"^{r]) contains two 
admissible arcs rjo and rji. Fix an index i. It follows from what is proved, that the inverse 
image S~^{vi) of any admissible arc rj, contains an admissible arc r/jQ C Hi^^. Applying the 
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same procedure to the curve rji^, we obtain an arc rjigi^ C S^"^{r]ig) f]Hi^. Finally, we get a 
curve Vmi-iru C ^-^^(i/ij f] 3-"'^+"' (Hi,) f]...f]Hi^. Then the statement of the lemma 
follows from the inclusion S^^rji^^i-^^^^i^) C Hi. 

Hence, for any point z & K there is a unique sequence 

i{z) = {. . . ,i_2,i-i,'io,H,i2, ■ ■ ■} 

such that S'^^{z) G Hi^ for any A; G Z. Due to Lemma [5] for any sequence i one may find 
the corresponding point z. Since the diffeomorphism S'^ is hyperbolic in a neighborhood of 
the set K the point z & K is uniquely defined by the sequence i{z). Therefore the set K is 
of the power continuum. The unit shift of the index to the left corresponds to the mapping 
5*^. The presence of this symbolic dynamics proves the theorem. □ 

9. Example 

Consider the equation x + px — qx = —a — 6sin(u;t) equivalent to the system 

y = qx — py — a — bsm{ut). 

We suppose that all the parameters of this system are real, and a, q and u are positive. 
Consider the vibro-impact system, defined by Eq. f|T3|) and the one-dimensional impact 
condition, where r = 1. The first component x{t) of the general solution of f|T3|) is of the 
form 

x = C+ exp(A+t) + C_ exp(A_t) + Lp{t). 

Here 



where 



_ -p+ ^/p^ + Aq ^ _-P- Vl^T^ 
A+ — > U, A_ — < U, 

(p{t) = — \- A sin(cjt) + B cos(wt), 
Q 

^ ^ biu"^ + q) ^ bpuj 



(w^ + g)^ + p^w^ ' {uj"^ + qy + p'^Lo'^ 

Note that ip(t) corresponds to the unique periodic solution of ffT3l) . The function ip(t) is 
positive for all t if and onlyif A^ + i?^ < a^/g^ or, equivalently, if < a^((a;^ + g)^+p^C(;^). 
The fundamental matrix of the corresponding homogeneous system, that turns to the unit 
matrix for t = is = 

1 /A+ exp(A_t) — A_ exp(A_|_t) exp(A+t) — exp(A_t) 
+ 4g \ -g(exp(A+t) - exp(A_t)) A+ exp(A+t) - A_ exp(A_t). 

The matrix A is equal to $(27r/a;). The spaces and are linear hulls of vectors 
col(l,A_) and col(l, A+) respectively and 

A+exp(A+t) — A_exp(A_t) 



JrC — A—, JrC — A-i_, JrC 



exp(A+t) — exp(A_t) 
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Conditions 4 and 5 hold true evidently. Clearly, ai2 > for all values of parameters of 
(fT3ll . Moreover, for all T we have > > i?*, that implies (jH]). So, the conditions of 
Theorem [3] are satisfied. 

10. "Soft" model and structural stability 

In this section we do not need the right hand side / of Eq. (|2]) to depend on the parameter 
H, so we omit this parameter in our notations. 

Suppose that Eq. (|2]) is defined for z E A. Fix a parameter value r G (0, 1] and denote 
Eq. dl]), corresponding to this fixed value of r by Note, that any Cauchy problem for 
the system with initial conditions z(to) = Zq, to G M, 

zq = col (xio, 1/10, ^o) e A 

has a solution. It is locally unique if xio ^ 0. In this case there is a value e > 0, such 
that on time intervals (to — ^,^0) and (to^^o + ^) the considered solution does not have any 
impacts and, consequently, coincides with a solution of Eq. ([2]). 

Consider a function X-{s) = 1 — H{s) where H{s) is the standard Heaviside step function, 
i.e. = 1 if -s < and X-{^) = if s > 0. For a fixed value r G (0,1] define 

a = — (logr)/7r. We define the function f{t,z) on the set M^"+^ \ A so that this function 
obtained is piecewise smooth on the whole space M^""*"^. Let > be a big parameter. 
Define 

/i(z/, z) = (— 2az/?/i — (1 + a'^)h'^Xi)x~ixi)ei. 

Here ei = col (1, 0, . . . , 0). The function h is piecewise continuous. 
In this section we study the following problems. 

1. When do the invariant sets of the shift mapping for Eq. (jl])^ persist provided the 
coefficients of Eq. ([2]) and/or the restitution coefficient r are slightly changed? 

2. When these invariant sets persist provided one replaces the impact condition with the 
perturbation hlu, z) corresponding to a big value of z/? 

Let g{t, z) = co\{gi(t, z), . . . , gn(t, z)) be a smooth mapping, T - periodic with respect 
to t. Suppose that it is small enough in the C° - norm together with dg/dz. Consider the 
following systems 

ik = Vk, Vk = fk{t, z) + 5(fc(t, z) and (14) 

ik = yk, Vk = fkit,z) + gk{t,z) + h{iy,z). (15) 

We denote the dynamical system, which consists of Eq. f|T4|) and Condition 1 where r is 
fixed by (JH])^^. Consider the solution Zrit,tQ,zo) of the system (jl])^ with initial conditions 
zito) = Zq and functions Zg^i,^rit,tQ, zq) and Zg^r{t,to, Zq), which are solutions of Eq. f lTSj) 
and (fT4|) _^ with the same initial conditions. Introduce shift mappings for Eq. f lT4|) _^ 
and ( !T5|) . by following formulae Gr{zo) = Zr{T,0, zq), Gg^r{zo) = Zg^riT,0, zq), Gg^u,r{zo) = 

Zg,ij,r(T, 0, Zq). 
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Theorem 6. Let vq G (0, 1]. Suppose that the mapping G = Gr^ has a hyperbolic invariant 
set K C (0,+oo) X M such that 0, 2:0) 7^ for all zq e K, t e [0,T]. Let U be a 

neighborhood of K such that U and the corresponding image G{U) do not intersect with the 
axis Oy. Let 

W = {t,z:te [0,T],z = z{t,0,Zo),Zo E U}. 
The following statements are true. 

1. For every e > there exists 5 > such that if r E {tq — 5, tq + 5) n(0, 1], 

dg 



max \q{t,z)\ < 5, max 

{t,z)(iW {t,z)&W 



< S, (16) 



the mapping Gg^r is well-defined in a neighborhood Uq of K . There exists a homeo- 
morphism hg^r '■ K — )■ Kg^r C Uq such that max \ hg^r{,x) — x\ < e, and the set Kg^r is 

hyperbolic invariant for the mapping Gg^r- Moreover, for any x E K 

hg^r{G{x)) = Gg^hgA^))- (17) 



2. Let Tq G (0, 1]. For any e > there exist Uq > 0, 6 > such that if v > Vq and 
conditions (fT6l) are satisfied, there exists a homeomorphism f]g^u.ro • K — )■ Kg^^ ro ^ 
such that max |?7„ — x\ < e and Kg^^ro o hyperbolic invariant set of Gg^y^^o- 

Moreover, r]g^,^^roiG{x)) = Gg^y^roiVg,i^,roi^)) f^f all x E K. 

Proof. Let us check item (1). The neighborhood U can be chosen so that any solution 
of Eq. ( fT4l) ,.^. which starts at the instant t = in the domain U, has at most a fixed value 
M of impacts. Let it be not true. Then there exists a sequence Zk E U, k E N of initial 
conditions, corresponding to solutions, which have at least k impacts over [0,T]. Without 
loss of generality, one may assume, that z^ ^ zq E U. If such a point exists for any choice 
of the neighborhood U, we may say that Zq E K. The solution 0, zq) has infinitely many 
impacts on [0, T], consequently, instants of these impacts have a limit point t* E [0, T]. Then 
z{t* , 0, Zq) = 0. This contradicts to assumptions of the theorem. 

The set U can be represented in the form U = Qo[J ■ ■ - [JQm, where every set Qk 
is compact or empty and consists of initial conditions corresponding to solutions of Eq. 
(fT^ ,.^. which have exactly m impacts on the segment [0,T], such that corresponding values 
of the normal velocity yi are always nonzero. Then the impact instants smoothly depend 
on the initial conditions. On every set Qk except Qq the mapping G is of the form G = 
IqoG^o. . .oG^^^oIi. Here /q maps the initial conditions 2:0 G A to the first impact instant t^, 
the point x{t^) E Hi and the velocity y{t^), corresponding to the first impact. The mappings 
G\ j = 1, . . . , k — 1, transfer the instant, position and velocity of the impact number j into 
correspondence to the same parameters of the j + 1-th impact. The mapping Ji transfers the 
instant, point and velocity of the last impact to the phase coordinates of the corresponding 
solution at the instant t = T. Clearly, mappings Jq, /i and G^ (j = 1, . . . , A; — 1) are G^ 
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smooth on the set Qi. Hence there is a a > such that the mapping G is smooth on 
the set U X [tq — a, tq + a]. 

The following auxiliary statement is analogous to the theorem on persistence of hyper- 
bolic invariant sets of diffeomorphisms [25]. The only difference is that we consider an 
embedding of a domain instead of a diffeomorphism of a manifold. The proof, given at [25] 
is still valid for the considered case. 

Lemma 7. Let d be a natural number, U dW^ be a domain. Let the smooth embedding 
(f) : U — 7- M'^ possess a compact hyperbolic invariant set K G U. Then for any e > there 
exists a 6 > such that if the mapping ip : U M.'^ is such that ~ 4'\\c^{u^R'i) < ^! then 
there is an embedding h : K satisfies the inequality \\h — id||f70(i^_^Rd) < e and 



for all X E K. Particularly, the set Ki = h{K) is hyperbolic invariant for the mapping ip if 
e is small enough. 

Applying Lemma [7] to the mapping Gr^ and its small perturbations Gg,r5 "we obtain 
that there exist 5 > and a neighborhood V d U oi the set K, such that for any r G 
(ro — 5, ro + 5) n(0, 1] and any perturbation g, satisfying conditions ffT6|) . there exists a 
homeomorphism hg^r-, topologically conjugating G and Gg^r in the sense of f lT7|) . Moreover 
hg^r =^ id as (5 — )■ and r — )■ tq. 

So, the mapping Gg^r is a diffeomorphism and the set Kg ^ = hg^r{K) is hyperbolic 
invariant. This proves the first part of the theorem. Let us start to prove the second one. 

Lemma 8. For ally^ > y^ > 0, any to G M and all functions f and g, satisfying conditions 
of TheoremlSi there exists a value vi > 0, such that if v > i^i, the solution of Eq. (|T5|l with 
initial conditions 

x{to) =xo e TTi, y{to - 0) = yo = col (?/i,o, Vo) : yi,o G {-y^, -y^) 

intersects tti at ti G {to, to + 2ti/v). Moreover there exists a A > 0, which does not depend 
on yo, such that for any v > Vi 



Proof. The transformation of the independent variable s = ut reduces (|T5|) for x ^ A 
to the form 



Here we denote by prime the derivative d/ds. The initial conditions (to,a^05?/o) are trans- 
formed to the following ones: 



h{(j){x)) = ij{h{x)) 



(18) 



yiiti) + royNfil < A/z/. 
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x'( + 2ax'i + (1 + a'^)xi — (/i(s/z/, x, vx') + gi{s/v, x, vx')))/v'^ = 0, 
x" — {f{s/h', X, vx') + g^sjv, X, vx'))/u'^ = 0. 



(20) 



so = uto, xi(so) = 0, x{so) = Xq, 



(21) 
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If V is big enough, Eq. ( l20l) on compact sets is a small perturbation of the system 

x'[ + 2cix\ + (1 + c?)xx = 0, x" = 0. (22) 

The first component of the general solution (p of Eq. ( !22l) is of the form 

Xi = Aexp{—as) sin(s + ip). 

The distance between neighbor zeroes of this function is always equal to vr, consequently, for 
u big enough and for any solution \1/ 7^ of Eq. (!20|) with initial conditions ( 12T]) the distance 
between neighbor zeros of the first component of \E' varies from tt /2 to 2it. For Eq. f|T5|) . we 
obtain that if the solution with initial conditions on U intersects with tti at the instant to 
and the corresponding normal component of velocity is negative. The solution intersects tti 
again at the instant ti G {to,to + 27r/z/). The ratio of normal components of derivatives of 
solutions of Eq. fl22l) corresponding to successive impact instants is —vq. Since the solutions 
between impacts continuously depend on initial conditions and parameters, there exists 
1^2 > 0, such that if z/ > 1/2, for any solution x{t) of Eq. (120|) with initial conditions ( !2T|l 
there is a A > 0, such that Ix'i^uti) + rox[{so)\ ^ A/z/^. This proves the validity of (fT9l) . □ 

Let the number m G {0, . . . , M} be such that Zq G Qm- If m = 0, the function 
Zg,u,rit,0, zq) is a solutiou of Eq. (ITSl) . Else, for u big enough and a. 6 > such that 
l^ol < the normal component Xi^g^,^^rit,0, zq) of the solution ^ ^,(^5 0, Zq) has exactly 
2m zeroes on the segment [0,T]. Denote them hy < 6^ < t'^ < 6'^ < . . . < t'' < 6'^ . 
The normal component ?/i,g,i/,r(^5 0, Zq) of the velocity of the solution Zg^,y^rit,0, zq) is posi- 
tive at the instants P and negative at the instants 6^. The values u and 6 may be cho- 
sen independently with respect to zq and m. The mapping G'g,i/,r is of the form Gg^^^r = 
Iq o o if^' o o if^' . . . o H^~^ o H^~^' o J]^. The mapping Jq transfers initial conditions 
Zq & K aX the time instant t = to the triple (t^,X^,y^). Here is the first zero of the 
normal component 0, 2:0)5 the tangential component = Xg^^^r{t^,0, Zq) and the 

value = yg^^^rit^ ,0, zq). Mappings H\ j = 1, . . . , k — 1 transfer triples (t\ X\Y^) of 
time instants P , impact points G tti and velocities , corresponding to the zero number 
2j — 1, to parameters 6\ X^' and Y^' , corresponding to the zero number 2j. The mappings 
H^' , j = 1, . . . , k — 1 transfer triples {6\ X^' , Y^') to {P^^, X^^^, Y^~^^) which correspond to 
the zero number 2j + 1. The mapping Ii transfers the triple [9^ ^X'^'^Y^') to the value of 
the corresponding solution for t = T. 

If u tends to infinity, 6 tends to zero, the instants P tend to the time instants of impacts 
of the solution of Eq. ([2]), corresponding to the same initial conditions, uniformly with 
respect to Zq, /o,i =^ /o,i in the metrics. Due to the result of Lemma [8] and Eq. (fT9|) . 
W W ^ id. It suffices to use Lemma [7] statement to finish the proof of Theorem 

El □ 

Remark. Theorem El applied to the results of Section 9 demonstrate that set of vibro- 
impact systems, satisfying Conditions 2 and 4 — 6, is non-empty and contains a subset open 
in X{J, n, T). 



REFERENCES 



21 



11. Conclusion. 

We have considered vibro-impact systems in their general form. It was shown that the 
existence of an unstable periodic motion which passes near the delimiter without having 
an impact may imply chaos. The corresponding sufficient conditions can be written down 
explicitly. This shows that grazing is not a single bifurcation, but a combination of two 
bifurcations that can coexist. This phenomenon takes place both for impulse and "soft" 
models of impact dynamics. 
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